Gapless excitations in strongly fluctuating superconducting wires 
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We study the low temperature tunneling density of states of thin wires where superconductivity is 
destroyed through quantum phase-slip proliferation. Although this regime is believed to behave as an 
insulator, we show that for a large temperature range this phase is characterized by a conductivity 
falling off at most linearly with temperature, and has a gapless excitation spectrum. This novel 
conducting phase results from electron-electron interaction induced pair breaking. Also, it may 
help clarify the low temperature metallic features found in films and wires whose bulk realization is 
superconducting. 
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Phase fluctuations of superconductors are responsible 
for a broad range of fascinating phenomena. Their effect 
is particularly dramatic in narrow wires, where the prolif- 
eration of phase slips induces a putative superconductor- 
insulator transition [H-Q . Experiments probing this tran- 
sition, however, challenge our understanding of the insu- 
lating phase, as they exhibit a low-temperature metallic 
behavior in thin wires where phase-slips proliferate 043 ■ 
Ref. Q discussed the phase-slip induced breakdown of 
superconductivity in a wire, speculating that a metallic 
phase arises. Could the strong phase fluctuations induce 
a finite quasi-particle (QP) density of states that main- 
tains a finite conductivity? 

Gapless quasiparticles are well known to exist in super- 
conductors when time reversal symmetry is broken [lCM — 
14| | . Gapless superconductivity, however, also appears 
due to proximity to metallic contacts where the order 
parameter is non- uniform 15|, [16[ . The pair breaking ef- 
fect occurs since the relative phase of the two electrons 
making up a Cooper pair (CP) gets randomized by the 
perturbation. This begs the question, whether fluctu- 
ating electromagnetic fields due to strong superconduct- 
ing phase fluctuations and electron-electron interactions, 
which arc manifestly non-uniform and introduce dephas- 
ing in normal systems |17| , can lead to the appearance of 
gapless superconductivity at finite temperatures as well. 

In this manuscript we study the low-temperature tun- 
neling density-of-states (tDOS) of the phase-slip prolif- 
erated regime. We argue that the scarcity of normal ex- 
citations and the blocking of the Cooper pair conduction 
channel give rise to strong dephasing, through electro- 
magnetic field fluctuations. This, in turn, should lead to 
pair breaking. From a self-consistent study of the tDOS, 
we find that at a broad temperature regime (roughly 
T < 0.1T C ), no hard spectral gap exists. Furthermore, 
because the conductivity is dominated by quasi-particles, 
it vanishes at most linearly in temperature, as opposed 



to an exponentially-suppressed conductivity characteris- 
tic of a gapped phase. Eventually, at very low tempera- 
tures, the QP are localizacd and the metallic phase ceases 
to be valid. These effects should be manifest in tunnel- 
ing measurements of wires with an increasing resistance 
upon cooling. 

Our argument follows from the dependence of dephas- 
ing on the dissipative response of diffusive electron sys- 
tems. For this purpose, it is insightful to interpret 
the response of the wires we consider in terms of co- 
existing normal quasi-particles and condensed Cooper 
pairs. In the phase-slip proliferated regime, where the 
normal resistance of a coherence-length segment obeys 
R.£ ~ Rq = h/4e 2 w 6.4fcf2, the conductivity is dom- 
inated by the normal QPs, as long as they remain dif- 
fusive. Similarly to normal diffusive systems, electron- 
electron interactions lead to the suppression of quantum 
interference of these diffusive QP after a typical dephas- 
ing time T^. Using the fluctuation-dissipation theorem 
the dephasing is dictated by the electrical response of 
the system, <r(T) [Hj]: 
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where on is the conductivity in the normal state, A is 
the cross-section area of the wire, and D is the diffusion 
constant. When the Cooper-pairs are formed but not 
condensed due to quantum phase slips proliferation, QP 
are expected to be scarce, and therefore <j(T) < on- This 
increases voltage fluctuations and the dephasing rate. 

An enhanced dephasing rate, however, may lead to 
the breaking of Cooper pairs, and therefore suppress the 
pairing gap. Indeed, if we assume a hard gap, A in 
the excitation spectrum we obtain a contradiction. If 
the QP density is exponentially suppressed n ~ e -A / T 
due to the pairing gap, then the dephasing rate is ex- 
ponentially enhanced. A strong dephasing mechanism 
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allows us to consider the effects of pairing on the tDOS 
perturbatively in the parameter Ar^(T) <C 1. Such a 
calculation yields a gapless excitation spectrum that ap- 
proaches the normal state tDOS (see Fig. [2] and detailed 
calculation below), in clear contradiction to the assump- 
tion of a finite excitation gap. This argument is valid 
at temperatures much lower than the mean-field T c yet 
above the localization limit of the QP. Note that thermal 
superconducting-phase fluctuations may also contribute 
to dephasing. These contributions are small, however, as 
QP's dominate transport in this regime. 
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FIG. 1: A qualitative phase diagram of a fluctuating wire, 
whose normal resistance is Rw, as a function of temperature. 
Here T c is the mean field transition temperature and Tq is 
the temperature at which quantum phase slips (QPS) dom- 
inate the response of the system. At t^t QPS proliferate, 
blocking the conductance of the superfluid channel. Gapless 
excitations appear at T < Tgap, where the resistance scales as 
a power law of the temperature, marked by an arrow. At Ti oc 
QP and CP are localized, and conductance is controlled by 
their thermally activated hoping with a typical gap E s . For a 
summary of all crossover temperature, and typical numerical 
values see Table I D For large the gap regime disappears. 
Here K is the Luttinger parameter, AF is defined in [2[| , and 
d is the Ginzburg Levanyuk number. The figure is not true 
to scale. 



Before we discuss the regime of interest, we first sum- 
marize the different temperature regimes of a strongly 
fluctuating superconducting wire, and their transport 
signatures. A qualitative phase diagram is depicted 
in Fig. [TJ Above T c [l9j], thermal pairing fluctuations 
and interference effects (Aslamasov Larkin and Maki- 
Tompson corrections), reduce the resistance. The for- 
mation of incoherent Cooper paris, however, tends to 
enhance the resistance. Together these reduce the re- 
sistance R{T) = R w — p-[(T — T c )/T c ] 3 / 2 , where R w is 
the wires's normal state resistance, and (T — T c )/T c <C 1. 
Near T c (|T - T c |/T c < Gi where Gi = (^^) 3/2 is 
the Ginzburg Levanyuk number) strong fluctuations con- 
trol the resistance. At temperatures below T c> a Cooper 
pair condensate forms, and shunts the normal excita- 
tions. Here thermally activated phase-slippage (TAPS) 



dominate, and the resistance follows an activation behav- 
ior, R LA mh(T) = ^g^e"^ 

At low temperatures T < Tq, where Tq is defined 
as A(Tq) = Tq, quantum phase slips (QPS) dominate 
the resistance of the superfluid. The behavior of super- 
conducting wires in this regime is dichotomized by the 
Luttinger parameter, K, which depends on see Eq. 
([5]) below. In wires with R^ < Rq, quantum phase slips 
are irrelevant and their resistance follows a power law 
temperature dependence. In this manuscript we focus 
on wires with R^ > Rq. At intermediate temperatures 
1/tkt < T < Tq, QPS are scarce but their density in- 
creases with reducing temperature, resulting in a power 
law temperature dependence, R(T) — R w (T / A) 2K ~ 3 . 
Here is the typical temperature at which phase 

slip proliferate, leading to a large phase slip fugacity: 
((1/tkt) = 1- At lower temperature , T < t^ t , QPS 
proliferate and the conductivity is dominated by diffusive 
quasi-particles. 

In the presence of a hard gap in the excitation spec- 
trum the number of normal quasi particles is exponen- 



tially small n c 



-A/T 



leading to an exponentially 



large resistance. However, the following remarkable cir- 
cumstances may lead to a metallic behavior character- 
ized by a power law resistance. The reduced density of 
normal excitation, resulting in the formation of Cooper 
pairs, increases the Nyquist thermal fluctuations of the 
potential. This finite temperature noise acts as a phase 
breaker for the superfluid, and consequently may lead to 
the vanishing of the gap in the excitation spectrum, a sit- 
uation known as gapless superconductivity. Unlike con- 
ventional superconductors, where the superfluid shunts 
the normal fluid, however, the proliferation of QPS block 
the superfluid channel and the resistance is dominated 
by these normal quasi-particles giving rise to a metal- 
lic behavior. The temperature at which the gap van- 
ishes, Tg a p, can be estimated from Ar^(T gap ) = 1. Us- 
ing Eq. ([TJ), for t$(T) and assuming the conductivity in 



the hard gap phase follows a(T)/a^ 



~ A / T , we find 



T gap «A/[ln(6oc/0 + lnln(6oc/0] 

The gapless regime, T < 
this manuscript. In this regime, as we outline below, the 
resistance of the wire follows a power law 



Tgap, is the main focus of 



a(T)/a N ~ (T/AY 



(2) 



Determining the power a requires summing the pertur- 
bation series in Ar^ to infinite order. Nevertheless, from 
Eq. © (cr(T) - t 3/2 T) together with the fact that 
T0 must diverge as T — > 0, we find that the conduc- 
tivity must follow a sub-linear temperature dependence 

a(T)/T > oo, corresponding to a < 1. We estimate 

a in two different ways: from the leading order correc- 
tion to the tDOS, and from the leading order correc- 
tion to the self energy. The former a ppr oximation gives 
<j{T)/<t n = v(T)/v Q w 1-(Atv»(T)) 2 |23|,[2i]. Substitut- 
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TABLE I: A list of all crossover temperatures, their defining 
relation, as well as their parametric and numerical values. For 
the numerical estimates we have used A ~ A(Tq) = Tq = 
0.9T C , b = 1, i?£ w 0.5-Rq. In addition, as the localization 
length satisfies R^ loc = 4_Rq, we have £ioc/£ = 8. We note 
that for these choice of parameters, T &s , p tkt > 1- 





Determined from 


Parametric 


Numeric 


Tq 


A(Tq) = T Q 




0.9 T c 


' KT 


C(tkt) = 1 


,-1 .. A 

T KT „njr_ Kc |-i/2 


0.2 T c 


Tgap 


AT</,(T gap ) = 1 


T ^p ~ ln iip +lnln «]p 


0.3 T c 


21 oc 




T loc « A (£/£ioc) 4 


0.0002 T c 



ing the expression for the dephasing time (p} and solving 
for the conductivity in the limit v(T)/vq <C 1, we ob- 
tain the power law a(T)/a^ « f^-^J- Alternatively, 
the leading order correction to the self energy, which is 
equivalent to a partial resummation of the infinite se- 
ries, gives rise to a sub-linear temperature dependence 
<j(T)/cjn ~ (T / A) . A discussion of this calculation 
and its validity appears at the end of the manuscript. 

The metallic behavior persists as long as the normal 
fluid remains diffusive. This breaks down at low tem- 
peratures, where the dephasing length exceeds the lo- 
calization length L<f,(T[ oc ) ~ £; oc . We estimate Xi oc by 
using the sub linear temperature dependence of the con- 



ductivity a(T)/(Ti\r 

\4 



(T/A) 2/5 in Eq. ©, and find 
A (£/£ioc)*) which for typical wires is well below 
T gap . Table I „ lists the different temperature regimes as 
well as their numerical values for typical experiments. 

The above argument shows that the assumption of a fi- 
nite dissipation mechanism in the low temperature phase 
of fluctuating superconducting wires holds self consis- 
tently. The finite conductivity derived from Eq. ([T]) is 
due to diffusive QPs, whose interaction dynamics gives 
rise to a fluctuating potential and hence dephasing. The 
dephasing rate not only gives rise to a finite conductiv- 
ity as described above, but also causes a finite tDOS. 
These effects can be probed in strongly fluctuating long 
superconducting wires whose resistance increases as the 
temperature is reduced, and whose total capacitance is 
large 
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261 ] . In this regime, the energy scale at which 
the tDOS reaches its maximum value is 1 /r^ (see Fig. [2]) , 
which according to these predictions should coincide with 
the value inferred from an independent measurement of 
the wire's conductivity following Eq. ((T|). 

Next, we prove the crucial point that if is indeed 
small, the QP spectrum cannot be gapped. We later 
use this calculation to estimate the emerging tDOS in 
the gapless regime. We carry out a calculation of the 
tDOS which is perturbative in At^. In the absence of a 
pair breaking mechanism, the perturbative correction di- 
verges, marking the opening of a pairing gap. Conversely, 
in the presence of strong dephasing, Ar^ <C 1, the QP ex- 



citation spectrum may be gapless. The tDOS is given by 
v t = -ilmG fl (r,r, e), where G R (r,r,e) is the retarded 
Green's function which can be expressed to second order 
in the pairing amplitude: G + G AG AGo(AA' 1 '). Here 
G (p,w Tl ) _1 =iw„-£p-H/(2T)Sign(a;„), G (p,o;„) _1 = 
iu n + ip+i/ (2r)Sign(cj T j) are the disorder averaged free 
Green's function in the vicinity of the second order phase 
transition, and A is the impurity ladder dressed ver- 
tex. We define Sv[e) = v{e) ~"° = -£lm I R (e), where 
I R (e) = I(iu> n — > e + iS) is the analytic continuation of 

I(u n ) = 27nsign(o; ra )r> |0 | ^ | ^ 9 | (3) 



q.Q 



(|2w n + n\+Dq 2 + T- L ) 



In order to describe correlations of the order param- 
eter in a superconducting wire we examine its micro- 
scopic action obtained from the BCS Hamiltonian by a 
Hubbard-Stratanovich transformation followed by an ex- 
pansion around the saddle point [27], HH . The low energy 
excitations of the system arc phase fluctuations whose ac- 
tion follow: 

S[<t>] = K/2 f dxdy {(d x cf>) 2 + {dyctf/N^} . (4) 



where y = v p t, N± = p^A/ir 2 and 
R = 4v AA 2 $ ^ 



2R{ 



(5) 



The partition function of the superconducting wire 
whose low energy excitations follow Eq. (|4]) has the same 
form as a classical partition function of an anisotropic 
two dimensional XY model. The system described by 
this model undergoes a Kosterlitz Thouless phase tran- 
sition between an ordered phase (superconductor) where 
quantum topological excitations known as phase slip arc 
bound in pairs and a disordered phase where phase slip 
pairs unbind [29| . Correlations of the order parameter in 
the disordered phase decay exponentially: 

(A(o:,T)A t (0,0)) = A 2 e- x/iKT e- T/TKT , (6) 

over a typical length ^kt, and time Tkt- This gives 



(AAt; 



q,n 



Aq^kttkt 



(i + g 2 £!- T )(i + ^iT) 



(7) 



Using Eqs. (J3|) and ([7]) we calculate the corrections 
to the tDOS in a fluctuating superconductor. The re- 
sult is shown in Fig. O We have assumed a prolifera- 
tion of phase slips which requires to be in the regime 
T < 1/tkt ^ A. The perturbative correction to the 
tDOS is small if A < l/r^, which in general is not sat- 
isfied in experiments. However, for illustrative purposes 
and to stay in the regime where our approximation is 
justified, Fig. [U shows the tDOS for Ar = 1,0.5,0.3, 
and the extension to the regime I/t^ < A is discussed 
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FIG. 2: The tDOS of a fluctuating superconductor given by 
Eq. ©. Here T/A = 0.1, and At K t = 2. Different curves 
are plotted for different dephasing time with At^ = 1, 0.5, 0.3 
corresponding to black dark gray and light gray, respectively. 
The tDOS approaches a constant value for short dephasing 
times. We note that the initial suppression of the tDOS at 
e > 1/ tkt is a result of the quasi-order of the superconducting 
phase at times shorter then the correlation time 1/tkt- 



below. In Fig. [2] one sees that the zero energy density of 
states decreases with increasing Ar^. As grows with 
decreasing temperature, we expect that v{e) develops a 
pscudogap, consistent with the power law temperature 
dependence of a. 

One might question the consistency of our calculation 
in the limit Ar^ <c 1, because the gapless tDOS we 
find can act as a shunt resistor, suppress quantum phase 
fluctuations, and restore local phase coherence [H-la, [3fl| . 
contrary to the assumed strongly fluctuating regime. 
This discrepancy is resolved since the phase-fluctuating 
regime (|7j) corresponds to a gapped kernel for the pair- 
ing field [31[. Hence, while gapless Fermionic excitations 
may introduce a dissipation term of the form R s | fi 1 1 A | 2 in 
the action for the pairing field, the substantial mass term 
in the action dictates the long time correlations render- 
ing the dissipation unimportant, and the system remains 
in the strongly fluctuating phase. Indeed, in the small f2 
limit, the pairing-field action (3lj with an R s > term co- 
incides with the Hertz-Millis action for the metallic phase 
of a strongly fluctuating superconductor [32 - 34 1 . In addi- 
tion, the phase diagram of dissipatively shunted Joseph- 
son junction chains also exhibits a disordered phase which 
is insensitive to small dissipation [IHs, Ho] ■ 

Let us now describe how to estimate the emerging 
tDOS and consequently the conductivity in the gapless 
regime, beyond the perturbative limit. Note that while 
the perturbative correction presented above is no longer 
small for a long dephasing time Ar^ > 1, it does not 
diverge and hence a calculation of the full density of 
states correction is possible. Similarly to the case of pure 
Coulomb interaction, is in principle due to higher or- 
der corrections in the perturbation series in the interac- 



tion strength. As such a calculation is beyond the scope 
of this letter, we use the dephasing rate obtained from the 
conductivity according to Eq. ^ in Eq. ©. Although 
it cannot be trusted, calculating the leading order to the 
self energy may give additional qualitative information 
about the tDOS in the regime Ar^ 3> 1, see details in the 
supplementray material. This leads to a finite sub-gap 
density of states at low but finite temperatures, which 
vanishes at T — > 0. In the limit Ttkt^Tt^ <C 1 we find 
a(T)/a a v(T)/l> q ~ l/(Ar (T)) = ^/^ loc T/A)^ 
were the last equality was obtained using ([T]) [23J, 12-41 ] . 
Similar models considering the role of a fluctuating pair 
correlator in a ballistic system in two dimensions revealed 
a QP peak in the spectral function [35|, [36| . 

In conclusion, we studied the tDOS of a supercon- 
ducting wire in the QPS proliferated regime. We found 
that in the disordered phase, associated with QPS prolif- 
eration, the conductivity has a sub-linear temperature 
dependence, and is dominated by quasi-particles with 
a substantial sub-gap excitation spectrum. This novel 
metallic phase may be related to the metallic behavior 
observed in low dimensional films and wires whose bulk 
realization is a superconductor. These predictions can 
be experimentally tested by tunneling measurements. In 
this regime, the energy scale at which the tDOS is ex- 
pected to reduce below its value in the normal phase, Uq 
is 1/t0, which according to these predictions should co- 
incide with the value inferred from an independent mea- 
surement of the wire's conductivity following Eq. ([T]). 
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Packard Foundation, BSF, DIP, the Alexander von Hum- 
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MICROSCOPIC PHASE ACTION 



In order to describe correlations of the order parameter in a superconducting wire we examine its microscopic action 
obtained from the BCS Hamiltonian by a Hubbard-Stratanovich transformation followed by an expansion around the 
saddle point [H, 0|. In the low temperature limit, this yields [H, 0: 

S = vqAA 2 , f L dx f 1/T dr { £ [In (p 2 ) - l] + 2£ 2 p 2 0' 2 + \tf + $ \p' 2 + ^(A \ , 
Jo Jo I * v 4, \ L V P J J 

where L and A are the wire's length and cross section, respectively, £q = 7r_D/8Ao, v p — \J (37r/2)Z?Ao the amplitude 
velocity, — y/irDAo {2AV c vo + 1) oc v p ^N± the phase velocity, V c the Fourier transform of the short range Coulomb 
interaction, N±_ = pj^A/n 2 is the number of one dimensional channels in the wire, vq the density of states, D the 
electronic diffusion constant in the normal state, and the SC order parameter is parameterized as A = Aope 1 ^, with 
Aq, the mean held solution. Rescaling the imaginary time by y — v p t, the low energy excitations of the system are 
phase fluctuations whose action follow: 



where the phase stiffness is 



S[<t>] = K/2 dxdy{(d x ^) 2 + (d v <j)) 2 /N ± } 



4isoAA 2 $ R Q 



K = 



2R{ 



(1) 



(2) 



The system described by this model undergoes a Kosterlitz Thouless phase transition between an ordered phase 
(superconductor) and a disordered phase where phase slip pairs unbind Q. Correlations of the order parameter in 
the disordered phase decay exponentially: 



(A(a;,T)A t (0,0)) = A 2 e- x/ ^ KT e - T ' TKT , 
over a typical length £,kt, and time tkt- This corresponds to 

A\^kttkt 



(AAt) 



(1 + ? 2 ^t)(1 + ^t) 



(3) 



(4) 



LEADING ORDER CORRECTION TO THE TUNNELING DENSITY OF STATES OF A 

FLUCTUATING SUPERCONDUCTOR 



The tDOS is given by 



v € = -~ImG R (r,r,e) = --Im f -^-G R (p,e), 

7T 7T J {2-Kr 



(5) 



where G R (r, r, e) is the retarded Green's function which can be expressed to second order in the pairing amplitude: 



2 



G{p,uj n ) = G (p,u n ) +T^2G {p,u n )A(q,uj n ,uj n + fi,)G (p + q,u n + Q)A(q,cv n + il, u n )G (p, w„)(AA t ) g! o- (6) 



Here: 



G (k + q,uj) 1 = i(ui) + — sign(o;) - £ 
a i ^ v 1 0(w(o; + fi)) 



2r|2cj + ft|+ZV + W 

and correlations of the order parameter are given by Eq. (|4]). The density of states is then given by 5v(e) — = 
— ilm I R (c), where I R {e) = I(ioj n —te + i6) is the analytic continuation of 

Z(w„) = 2ot signK)T^— _! ■ (7) 

^(|2w„ + O|+Dq 2 + r ) 2 

Using Eq. (|4]) to describe the phase fluctuations in a phase-slip proliferated wire, in the low energy limit -C t^t 
we may approximate Eq. (J2J) as 



/(w„) w 27risign(w n )— — — — zjrrT^ n2 2 / 

(|2w„| +r/) 2 V l + n 2 7^ T 7 

7risign(w„)Ag f i 



where ^(z) is the digamma function 



2 2ttT 2ttTt kt J \2 2ttT 2ttTt kt 



LEADING ORDER CORRECTION TO THE SELF ENERGY 

The leading order correction to the self energy is given by G _1 = Gq 1 — S with: 

E = E T E 6 ( fc + «' w + fi)(AA t ) 9 , n A 2 (w,w + n,g). (9) 
q a 

The integral over fermionic momentum is dominated by £ w 1/r. Since, ujt,Q,t, Dq 2 r <C 1, we can approximate 
G(fc + q>, w + fi) w G(fc, w). This gives 

E « G(^)VrV ,„ 9 ("(" + ")) g <AAt), n 
« n 4r + fi l + + 1 / T ^ 
= G(k,u)A{w). (10) 

Using this expression for the self energy we can write the green's function as: 

Gik.uj)- 1 = i(w) + ^-sigiuV) - & - E(w) (11) 
= iw-f-T^rrACw), (12) 
where cj = u; + ^rsign(o;). The density of states is given by: 

v(ioS) = — J dkG(k,Lj) = —v$ j d£— 



= * /A s (13) 

^w 2 + ^(w) 
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where the odd integral over £ vanishes. In the limit of lot <C 1 we have: 



v{ilo) = Vq 



isign(cj) 



Vq- 



«sign(w) 



y/1 + 4T 2 A(CJ) y/AT 2 A(i0) 



(14) 



where the last approximation is valid beyond the perturbative limit where 4t 2 A(lo) 3> 1. 

In order to evaluate 4t 2 A(lo) : we note that At 2 A{lo) — 2 ^si°-n(ij) wnere I(lo) is given by Eq. 0. Using Eq. 
the limit Ttkt *C 1 we find: 



4t 2 A{lo) 



A 2 



2 (|2w| + l/r^)2 l4vr 
A 2 



1 



2 2ttT 2itTt kt ' 



1 w 



2 2vrT 2ttTt kt ' 



coth( 



2Tt 



-) 



1 



- — [«7r + 2it:Ttkt — 2wtkt] + 1/2 



2 (|2w| + l/r ) 2 Utt 

Here we have assumed to ~ T <C 1/trt. Performing the analytic continuation i<j — > e + i<5 we find 

1 



At 2 A{ilo -> e + id) 



A 2 



[*7T + 2mTr KT - 2eT KT ] + 1/2 



2 (-2ie + l/r ) 2 ^ 4tt 
T 1 (-2 l£ + l/r,) 2 (i _ i £TKT _ ^ TA ' T 



The density of states is given by 

v(e) — Qv(iu) — > e + iJ) = 3 



(-2ze + l/r ) 



Jl/2 - ±€T KT - Tt kt 



In the low temperature limit Ttkt, Tt^ -C 1, we can replace v(T) = — J dev(e)-J^ ps i/(e = 0, T), leading to: 

i/(T) _ 2^2 
~W ~ Aor^(T)' 



(15) 



(16) 



(17) 



(18) 
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